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General Relativity's Kerr metric is famous for its many symmetries which are responsible for the 
separability of the Hamilton-Jacobi equation governing the geodesic motion and of the Teukolsky 
equation for wave dynamics. We show that there is a unique stationary and axisymmetric Newtonian 
gravitational potential that has exactly the same dual property of separable point-particle and wave 
motion equations. This 'Kerr metric analogue' of Newtonian gravity is none other than Euler's 18th 
f***) , century problem of two-fixed gravitating centers. 
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INTRODUCTION 



The Kerr metric describing rotating black holes, apart from its extreme physical significance is also abundantly 
enriched with symmetries [l|. The existence of these symmetries is mirrored in the equations governing the motion 
of particles and waves in the Kerr gravitational field. The Hamilton-Jacobi equation for particle geodesic motion 
and the Teukolsky equation for wave dynamics are known to be separable, thus offering a deep understanding of 
the underlying physics Q. Indeed, the impressive theoretical advance in black hole physics that took place in the 



black-hole 'golden age' 1970s-1990s was based on this special property of the Kerr metric jl| 



On the other hand, Newtonian gravity of course cannot create black holes. However, some Newtonian gravitational 
potentials may possess many of the symmetries of the Kerr metric. For instance, the notion of a 'third integral' of 
motion (besides the obvious integrals of energy and angular momentum along the symmetry axis) is a well-known 
concept in Newtonian gravity, e.g., in the context of galactic potentials Q, and is associated with integrability, that 
is, the particle motion having as many conserved quantities in involution with each other as the number of degrees of 
freedom In black hole physics this integral of motion is known as the Carter constant [j| and its existence is linked 
' with the 'hidden' symmetry of a rank-two Killing tensor in the Kerr metric (in contrast with the more obvious 
symmetries of stationarity and axisymmetry which lead to a conserved orbital energy and angular momentum along 
* ■ the spin axis, respectively). 

A famous example of a Newtonian potential endowed with a Carter-like constant is the Euler solution for the 
£f) • two-centre problem, i.e. the gravitational field sourced by a pair of spatially-fixed point masses. A recent account of 
t-H ' the integrability of particle motion in the Euler field and a simple derivation of its Carter constant has been provided 
by Lynden-Bell Q. More recently, Will [§} revealed an unexpected aspect of Euler's 250-year old solution. Within 
the whole class of solutions of the Poisson equation, there is a unique axisymmetric and reflection-symmetric solution 
(besides the spherically-symmetric monopole solution) which has a third integral like the Carter constant. This special 
potential is none other than the Euler solution (with equal- mass centres). Remarkably, the mass multipole moments 
in which the Euler potential can be decomposed are interrelated in exactly the same way as the mass multipoles of 
the Kerr metric! 

Based on these results one is compelled to argue that the Euler potential is the Kerr analogue in Newtonian gravity 
(this analogy was first noticed long time ago, see Refs. [9j,[l0|). However, the close relationship between the two 
gravitational fields has been established within the context of point-particle motion, i.e. the presence of a Carter-like 
constant and its relation to the separability of the Hamilton-Jacobi equation [ll[ . Knowing that the Kerr metric is 
also special with respect to wave dynamics [2|, the emerging question is whether we can push the Euler/Kerr similarity 
even further, by comparing these two gravitational fields with respect to wave-dynamics properties. 

More specifically, a key question to be addresed is this: does the Euler problem lead to a separable wave equation 
in the context of Newtonian gravity, in the same way that the Teukolsky equation for the Kerr metric has a separable 
solution? Moreover, is this the only Newtonian potential (besides the spherically symmetric one) with the dual 
property of rendering separable both the Hamilton-Jacobi and the wave equations? If the answer is affirmative, then 
the analogy between the relativistic Kerr metric and the Newtonian Euler potential would be 'complete'. Besides 
having significant theoretical value in itself, such an analogy between the two fields could be used to get a deeper 
insight into these two fields or even, perhaps, predict new properties for them. 

In this paper we provide an affirmative answer to the above questions. That is, we show that the Euler potential 
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is the unique axisymmetric and reflection-symmetric Newtonian potential that allows the separability of the scalar 
wave equation (in addition to that of the Hamilton- Jacobi equation), using a spheroidal coordinate system. 

The rest of the paper is organized as following: In Section [2] we discuss the occurrence of the Carter constant in 
Newtonian gravity and the connection between the Euler potential and the Kerr metric. In Section [3] we formulate 
the wave equation in Newtonian gravity and discuss its separability. In this section we show, with a modicum of 
effort, that in usual spherical coordinates the separability of the wave equation is possible only for a monopole (i.e. 
spherically symmetric) potential. Within the more general class of prolate/oblate spheroidal coordinates though, the 
wave equation is separable if and only if the potential is the Euler one. Finally, in Section|3]we discuss the similarities 
between the basis functions of Poisson's equation for the Euler problem and the corresponding ones for the wave 
equation in the Kerr metric. 



2. EULER'S 'TWO-CENTRE' GRAVITATIONAL FIELD 

2.1. The analogue of Carter's constant in Newtonian gravity 

Particle motion in stationary and axisymmetric gravitational fields is characterized by a conserved energy E and 
angular momentum L z along the symmetry axis. In this case the connection between the symmetries of the background 
field and the nature of the conserved quantities is, more or less, obvious. The existence of an additional third integral 
is a more rare property, and is not related to any obvious symmetry of the gravitational field. In General Relativity, 
and more specifically in the case of the Kerr metric, the presence of such a conserved quantity-the Carter constant-is 
a consequence of the underlying presence of a rank-two Killing tensor. The Carter constant itself is known to be 
the sum of the total angular momentum L 2 plus a spin-dependent term; its true meaning (to some extent) has been 
clarified only recently [12| . 

In Newtonian gravity there are very few potentials V that are solutions of the vacuum Poisson equation W 2 V = 
and which are endowed with a third integral. Within the general class of stationary and axisymmetric potentials 
the simplest one possessing a third inegral is the superposition of a monopole-dipole field. Using standard spherical 
coordinates {r, 9, ip}, this potential is (we set G — c — 1 adopting the geometrized units of General Relativity) 

M d • r 

Kld(r) = __ + _ (1) 

where M is the mass of the source and d is its dipole moment. According to the Landau & Lifshitz Mechanics 
textbook [11] . V^nd is the most general potential in spherical coordinates that leads to a separable Hamilton- Jacobi 
equation for particle motion 1 . 

Another example of a potential with a third integral constant is provided by the famous solution found by Leonhard 
Euler in 1760 for the problem of two fixed gravitating masses. This system consists of two static point masses mi 
and m 2 with a fixed distance 2a between them. Without loss of generality we can assume that the two masses are 
located on the z-axis (which thus becomes the symmetry axis) at symmetrical places around the origin. The resulting 
potential is 

t r mi m 2 t r mi m 2 , . 

V E = -- -- -> V E = (2) 

|r — az\ |r + az\ r% r 2 

where ±az are the positional vectors of the two masses with respect to the origin (a 'hat' denotes a unit vector) and 
ri,r 2 are the distances of the given observation point from the two masses. 

The third integral associated with point-particle motion in this potential is then 

Q = -hi • h 2 + (m 2 f 2 - mifi) ■ za (3) 

where f j is the unit vector from the i-th mass to the particle and hj is the particle's reduced angular momentum with 
respect to the i-th mass; i.e. h; =r;xv with v the particle's velocity. It should be noted that the explicit form of 
Q can be derived in more than one ways: for example, Lynden-Bell [7| provides a straightforward derivation of this 
result, while Ref. (Tl| arrives at the same result by separating the Hamilton- Jacobi equation. 



1 Interestingly, on page 901 of Ref. |l4j ] it is stated that the integrability of motion in the potential V m d, inspired Misner to suggest Carter 
to look for a third integral in the Kerr metric. 
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In a recent paper [8], Will was able show that there is a unique 2 stationary /axisymmetric and reflection-symmetric 
Newtonian potential that possesses a third integral that is quadratic with respect to the momentum (as the Carter 
constant). It turns out that this special potential is identical to the reflection-symmetric case of the Euler solution 
@, that is when mi = m 2 = M/2. For later convenience this potential will be denoted as V s e- Hence, 

_ M/2 M/2 

Note that the additional requirement of reflection symmetry has eliminated the potential V m d (eqn. (JlJ) from the list 
of potentials with a third integral since the dipole field lacks this property. 

Will arrived at the potential V s e by asking a question with bearing on the Kerr metric. Considering a general 
decomposition of a stationary, axisymmetric, and reflection-symmetric potential in multipole moments Mi (where 
I = 0, 1, 2, ... is an integer), 

V = -fl^T P ^9), (5) 

what is the relation (if any) obeyed by M 2 /, in order for V to lead to particle motion with a third integral/Carter-like 
constant? It turns out that a recurrence relation does exist between the moments and is given by 

M? ( Mo \ 1 

M t = — M1-2 -> M 2 i = M [ — - ) , (6) 

where Al = M, M 2 is the quadrupole moment and so on (only the even-order moments M 2 ; appear here, while the 
odd-order moments are all vanishing, A/ 2 z+i = 0, as a result of the imposed reflection symmetry). 

Remarkably, this is the very same relation obeyed by the mass multipole moments of the Kerr metric, as defined 
by Hansen and Geroch [l5j . Exactly this relation, together with a similar one for the current-mass multipoles (which 
of course have no analogue in Newtonian gravity), completely characterizes the Kerr metric, and encapsulates the 
mathematical form of the famous 'no- hair' theorem for Kerr black holes [l4j |. 

The nature of the Newtonian source with the structure of eqn. © is easily revealed if we use in (JSJ) the following 
identity for the Legendre polynomials Pi (/i) (correcting a typo in [8j , in the generating Legendre function formula just 
above eqn. (22) of that paper) 

00 

E t2lp M = 



2 

Then ([5]) reduces to the potential V s e of (0| with 3 




-1/2 



(7) 



For a prolate gravitational field, M2 > 0, the distance a is a real number and the resulting gravitational field source 
can be identified with Euler's two-centre system (two point masses mi = m 2 = M/2 located at z = ±a). On the 
other hand, the oblate field M 2 < (which is the one that corresponds to the multipolar structure of the Kerr metric) 
entails an imaginary distance 4 ia. The mass moments then become 

M21 = Ma 21 , (prolate) (9) 

M 2 i = M(-a 2 ) 1 , (oblate). (10) 



2 The uniqueness of this potential has been formally proven very recently, see Ref. [l3t | . 

3 Based on intuition alone, one might have expected that the Newtonian potential closer to the structure of the Kerr metric would be the 
field of a circular ring distribution of matter, as a substitute of the ring singularity of Kerr black holes. However, the explicit calculation 
of the multipoles M21 = ( — 1) Mb 21 (21 — 1)!!/(2 ; Z!) (where b is the ring's radius) reveals a clearly non-Kerr multipolar structure. This 
result is probably due to the fact that the usual spherical coordinates used to decompose a Newtonian potential into multipole moments 
are not exactly analogous to the Boyer-Lindquist coordinates in which the Kerr metric is usually described. 

4 In order to avoid any misunderstanding we will keep the distance a real, and add the imaginary unit i whenever we discuss the oblate 
case. 
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The connection between the Euler potential and the Kerr metric is clearly intriguing although there are still some 
distinct differences. The Kerr metric describes an oblate gravitational held, as one would naturally expect due to its 
intrinsic spin. In contrast, the real space Euler potential is prolate by construction. This difference is reflected in the 
sign of the quadrupole moment M 2 ; Kerr's quadrupole moment being negative, while Euler's quadrupole moment is 
positive. 

Opposite sign quadrupole moments would lead to a qualitatively different particle motion in the two fields. This 
difference can be annulled artificially by 'rotating' the distance a on the complex plane, i.e. a — > ia. Although 
this procedure cannot be produced by a suitable rearrangement of the masses in the physical space, the resulting 
potential is a legitimate, real- valued solution of the Poisson equation. This fact is obvious from the expansion ([5]), 
but it can also be derived by considering the symmetric Euler potential Q and changing a to ia. Thus, after all, it 
is more appropriate to think of the oblate Euler potential (which has the same set of mass moments as Kerr) as the 
'Kerr-surrogate' in Newtonian gravity (a more detailed discussion on the extraordinary similarities between the two 
fields will be presented elsewhere [l7|). 

The oblate Euler solution can be interpreted as a two-centre system in the complex a space with the masses located 
at ±iaz. We should note that the transformati on a —t ia should be accompanied by expressing vector amplitudes 
(denoted as | . . . |) according to the rule |k| = Vk ■ k. This product is a complex function and in order to keep the 
square root single- valued square root we should place a branch cut along the negative real axis of the vector inner 
product. 

With this complex-plane extension of the Euler solution, the general potential Ve assumes real values only in the 
symmetric case m\ = m 2 since then the two terms in @ are complex conjugate to each other, |r — iai\ = |(r + iaz)\*. 
After some lines of algebra the resulting oblate potential, henceforth denoted as V s e, takes the real- valued form 

V sE = -^VR 2 +r 2 -a 2 , (11) 

where R 2 = \J [r 2 — a 2 ) 2 + (2ar • z) 2 . 

A similar result holds for the Carter-like constant Q, eqn. ([3]): assuming an oblate Euler field, the vectors r 2 and 
f 1 are complex conjugate with respect to each other; this makes f 2 — f 1 a purely imaginary quantity. The same is 
true for the inner product hi • h.2, since we have that hi = h^. Given these properties, and despite the presence of 
the imaginary parameter ia, the constant Q is a real number provided the system is symmetric, m\ = m 2 . 



2.2. The separabillity of the Hamilton-Jacobi equation 

Before moving to the main topic of this paper (the separability of the Newtonian wave equation) it is worth 
considering the separability of the Hamilton-Jacobi equation for point particle motion. The aim here is not to provide 
any new results but merely highlight some key points with relevance to our discussion. 

As discussed in more detail in Ref. [ll|, the separability of the Hamilton-Jacobi equation in the gravitational field 
of the Euler potential can be achieved in spheroidal coordinates {£, rj, ip}. The following discussion can be made more 
compact if we simultaneously study both the prolate and oblate Euler fields. To this end we will have to use on the 
same footing both prolate and oblate spheroidal coordinates: all formulae will be expressed in both coordinates; the 
upper (lower) sign will be reserved for prolate (oblate) coordinates, whenever such a discrimination is needed. 

The spheroidal coordinates {£, rj, ip} are related to the standard spherical {r, 9, ip} and cylindrical {zu, z, ip} coordi- 
nates as 

^=* 2 Tl)(l-') 2 )- z = r tri, V = V (12) 

and 

.-=,„p ± (,*-i)r, ^= ^ ± j"_ 1)r («) 

where ro is a constant length-scale. The spheroidal coordinates span the domain 1 < £ < 00 (prolate coordinates), 
< £ < 00 (oblate coordinates), — 1 < rj < 1 and < ip < 2tt. In other words, £ can be viewed as a radial coordinate, 
rj as the cosine of a latitude, and <p as the usual azimuthal angle. 

The spheroidal system at hand is in a sense 'adapted' to the two-centre problem. In terms of the radial distances 
r*i , r 2 from the two centres and identifying ro = a we have 

£ = — , r, = — prolate 

la la 



r\ +r 2 r 2 - T\ 



2a ' 2ia 



(oblate). (14) 
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Relating the spheroidal coordinates to the spherical coordinates used in the earlier discussion of the two-centre system 
is a straightforward task in the prolate case. The oblate case, on the other hand, is somewhat more involved as a 
result of using the imaginary distance ia. The coordinate relation is 



R 2 + (r 2 - a 2 ) 



2a 2 



r] = sgn(cos 0) 



R 2 - {r 2 



2a 2 



(15) 



where sgn(x) is the usual sign function and R was defined immediately below eqn. (|11[) . It is noteworthy the fact that 
the oblate system contains points where £ = 0. Given that r\ — r?>, we can have r • z = and r 2 — a 2 < 0. Hence on 
the whole equatorial disk region r < a, r\ and r 2 are purely imaginary complex conjugates and £ = 0. 

According to Ref. [TTJ] the most general stationary and axisymmetric potential V({;, if) that renders separable the 
Hamilton- Jacobi equation in spheroidal coordinates (referred as "elliptic coordinates" in the textbook) is 



^scp 



£ 2 =F if 



(16) 



where F, G are arbitrary functions. Although Ref. fTTjj obtains this result in prolate spheroidal coordinates, it is 
actually true in oblate spheroidal coordinates as well. In the ensuing analysis of the Hamilton- Jacobi equation the 
Carter-like constant Q emerges as one of the separation constants. This result is consistent with the third integral 
derivation of Ref. Q • 

If, in addition, the potential V scp is required (as it should!) to solve the Poisson equation we then arrive at the less 
general expression (see Appendix [A] for the derivation) 



Vscp — 



M + Bi] 
S, 2 T V 2 



(17) 



where A,B are constants. With the help of (|14|) we can see that, in fact, this potential is identical to the general 
Euler potential, given by eqn. @, with the identification 



.4 
.4 



mi + m 2 
a 

mi + m 2 



B 
B 



mi — m 2 
a 

. mi — m-2 



(prolate) , 
(oblate). 



(18) 



Once again we can see that only the symmetric system mi = m 2 generates a real- valued oblate potential (B = 0). 
This potential then automatically becomes reflection symmetric with respect to the z = plane. 



3. THE WAVE EQUATION IN NEWTONIAN GRAVITY 



3.1. Formulation and separability 



As we emphasized in the Introduction, the Kerr metric is special because (among other things) it allows the sepa- 
rability of both the Hamilton- Jacobi and Teukolsky partial differential equations, in the well-known Boyer-Lindquist 



coordinates [14{. Our objective of this Section is to show that the Euler gravitational field is equally special, in the 
sense that it is the unique Newtonian potential (with the symmetries of stationarity, axisymmetry, and reflection- 
symmetry) that leads to separable Hamilton- Jacobi and wave equations. 

Before showing this in detail, we first ought to clarify what we mean by 'Newtonian wave equation'. We define this 
equation as 5 

= -KFf (19) 



5 An alternative definition of the wave equation could be based on the (time-independent) Schrodinger equation - it can be shown that 
also this equation is separable in the Euler field. However, we prefer working with our definition of the wave equation as it is second 
order with respect to the time-derivative and it does not contain the Planck constant, which would otherwise be needed to balance the 
dimensions in the corresponding Schrodinger equation. We thank the anonymous referee for pointing out this possibility. 
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where ty(t, r) is a (massless) test-scalar field which we assume to be artificially coupled to the gravitational potential 
V and k is a suitable constant so that the equation is dimensionally balanced. Assuming hereafter a harmonic 
time-dependence 'J ~ e lujt we have the corresponding time-independent equation 



V * + I 72" _ kV 1 * = °- 



(20) 



For the moment we will use the usual SI units in order to clarify some yet undefined issues with respect to k, but later 
on we will return to the geometrized units which we employed previously in the paper. Since a Newtonian gravitational 
potential has the dimensions of velocity squared there are only two possible choices for k: (i) k = c 2 /(GM) 2 and (ii) 
k = w 2 /c 4 , plus combination of those two like (c/GM)(uj/c 2 ). If we introduce an extra length scale in the problem, 
like later on when we will study the separablity of the wave equation in spheroidal coordinates, another possible choice 
for k will show up; k = l/(a 2 c 2 ). Furthermore, combinations of the latter expression for k with the previous ones, 



like k = uj/(ac 3 ), could also correctly balance the dimensions in eqn. ([20)) . However, one has to bear in mind that all 
a-dependent choices are pathological in the limit a — > 0, because then the potential term is divergent and we end up 
with a singular wave equation in a Keplerian potential. A similar situation arises for the choice k = c 2 / (GM) 2 ; the 
resulting wave equation is problematic in the limit of a freely propagating scalar field, i.e. M — > 0. 

Hence, the only viable, singularity-free, choice is n — uj 2 /c 4 and is the one we will adopt in the following analysis 
despite the fact that it alters our picture of a fixed background gravitational potential. At this point, after having 
clarified the various possibilities for k, we can resume the use of geometrized units c = G = 1 (and set k = to 2 when, 
later on, we compare the Newtonian wave equation in spheroidal coordinates with the corresponding equation for 
Kerr). 

We begin our analysis of eqn. (|20[) by first attempting separation of the wave equation in spherical coordinates. 
Using the multipolar expansion ([5]) it is easy to see that only the monopole term —M/r allows separabality. Hence 
from early on we are forced to consider the problem in alternative coordinate systems. 

Motivated by the separability of the Hamilton- Jacobi equation in spheroidal coordinates, it makes sense to study 
the wave equation in the prolate/oblate coordinate system {£, r), ip} (having set tq — a as before). In these coordinates 
the Laplacian takes the form (recall that the upper (lower) sign corresponds to prolate (oblate) coordinates): 



V 2 * 



2(£2 zp ^2) 



dd(Z 2 t + d v {(i -v 2 )d^} 



d 2 v 



aTTl)(l-') 2 )' 



(21) 



We then look for a separable solution of the form 



= R{()S{T])e lmip . 



(22) 



Inserting this in ([20|) we obtain, 



_ 2 d 2 R dtR m 2 
+a 2 K 2 TT)V-^) = 0. 



(23) 



Clearly, as far as separability is concerned, the key term is (£ 2 =F rj 2 )V. The most general functional form of V that 
accomplishes separability is obviously 



(eTV 2 )V = F(0 + G( V ), 



(24) 



where F, G are arbitrary functions. This is nothing else but the potential V sep of eqn. (fl"6)) that separates the Hamilton- 
Jacobi equation as well. We then end up with the following pair of ODEs: 



9„ {(1 - r? 2 )V?} - 
(£ 2 T l)fl£iZ + 2^i? 



1 — t]' 



± (ujari) 2 + Ka 2 G{ii) + E 



5 = 



na 2 F(0 - E 



R = 



(25) 
(26) 



where E is the angular separation constant. 
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Now, if we impose the additional requirement V 2 V — 0, so that V corresponds to a gravitational potential, we find 
(see Appendix \K§ that separation of the wave equation is achieved only for the Euler gravitational potential (and in 
the case of oblate coordinates only for the symmetrical one, if we insist on real- valued potentials): 



(mi + m 2 )i + (m 2 - mi)rj 

V = V E = 775 jt prolate) 

a (s — V ) 

V = V sE = - 2 2^ (° blate ) > ( 2? ) 



where M = 2mi = 2m 2 in the oblate case. 

In combination with the discussion of Section [2j we have thus shown that the Euler field Ve (and V s e), besides 
being the unique stationary and axisymmetric potential with a Carter-like constant and with a Kerr-like multipolar 
structure, it is also the only one that leads to fully separable Hamilton- Jacobi and scalar wave equations in spheroidal 
coordinates. Based on these special properties we are compelled to conclude that the Euler field indeed qualifies as 
the Kerr analogue of Newtonian gravity. 

It is also interesting to determine the multipolar structure of Ve (or V s e) in the spheroidal coordinate system 
(we have already discussed the multipoles in spherical coordinates, see Eq. ©). To accomplish that we start by 
considering an axisymmetric potential of the form 

V(t,rj) =XX(0<Sn(»7) (28) 

n 

Then from 'V 2 V = 0, written in spheroidal coordinates (either prolate or oblate), we obtain a pair of Legendre-type 
equations 

(1 - rf)d 2 S n - 2ry9„S„ + X n S n = 0, (29) 

and 

(C 2 T l)dpl n + 2Cd s TZ n - \ n K n = 0, (30) 

where A n = n(n + 1) is the corresponding angular eigenvalue. The solution that is reflection symmetric and asymp- 
totically vanishing at infinity is (after a trivial rescaling of the summation integer) 

oo 

/J Q!2n Q2n (Q f*2n (v) (prolate potential) 
V = <! (31) 

oo 

at2 n iQ2n(iQP2n{'n) (oblate potential) 

n=0 

where are constant coefficients and Pk,Qk are Legendre polynomials of the first and second kind, respectively 
[l6|. The extra i appearing in the oblate potential has been added to ensure a real- valued expression, assuming the 
coefficients ai n are real (recall that Q 2n is purely imaginary if its argument is purely imaginary). Eqn. (|3ip is the 
analogue of the multipolar expansion ([5]) but in spheroidal coordinates. The correspondence can be formalized if we 
use the large argument approximation for the Legendre polynomials either for real or imaginary argument [Trl ] 

Q^»)»t^r™. (32) 

It then follows by a suitable dimensional rescalling of the terms in the series (£ is dimensionless) that 

f -Er=o^2nK)" (2 " +1) - p 2nW (prolate potential), 
V(Z^oo,t,)m { (33) 
I -EZ a M 2 rM)- {2n+1) P2n(v) (oblate potential), 

where the spheroidal multipole moments arc then defined as 

(2n)!a 2 " +1 , x 

M ^ = ~-^TW^. (34) 
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As we show in Appendix [B] both prolate and oblate series share the same coefficients «2n , namely, 

a 2 „ = -(4n + l)— . (35 
a 

Hence, the spheroidal moments of the Euler potential are given in closed form by 

(2nV 

M 2n = ,. K u ,, Ma 2n , n = 0,1,2,... (36) 
(4n — 1)!! 

with M.q = M, A^2 = 2Ma 2 /3, etcetera. Most noteworthy, these moments do not obey the Kerr rule © and they are 
all positive regardless of prolateness/oblateness. The latter property is not too surprising given that the two types of 
potential are build on a different set of base functions (<?2z(£) or iQ2i(*£); respectively), which means that the same 
coefficient terms do not mean identical potentials. At the same time the deviation from the Kerr multipole moments is 
due to the fact that the Newtonian multipolar expansion is coordinate-sensitive, in contrast to the coordinate-invariant 
formulation of the Geroch-Hansen relativistic moments. 



3.2. Wave equation: Euler vs Kerr 



Having established the separability of the wave equation we now can focus on a closer comparison of that equation 
against the Kerr scalar wave equation. Based on our previous discussion, it makes sense to consider for this comparison 
the oblate reflection-symmetric Euler potential V s e which is the one closest to Kerr as we argued above. 

With the identification V = T4e, the angular and radial equations ([25]) and ([26]) become 



d n {(l-r, 2 )d n S} 



1 — r\' 



E 



(£ 2 + l)df i? + 2£dsR 



(aw£) 2 + 



5 = 0, 



kMo£ - E 



R = 0. 



(37) 



(38) 



The first equation has exactly the same form as the angular equation of the spin s — (scalar) spheroidal harmonic 
functions - an equation produced during the separation of the scalar wave equation (i.e. the s = Teukolsky equation, 
see eqns. (4.9) and (4.10) of 0) in Boyer-Lindquist coordinates {t, r B L, 9bl, <Pbl}- Matching the two equations requires 
that we identify r\ = cos#bl, where #bl is the Boyer-Lindquist latitude coordinate. 

The second equation resembles the corresponding s = radial Teukolsky equation without, however, being identical 
to it. We can nevertheless push the resemblance further by replacing the dimensionless coordinate £ in ([38)) with the 
following function of the Boyer-Lindquist radial coordinate tbl^ 



£ = tan$(r BL ), 



with 



$(r B L) 



dx 



A(a;) = x 2 - 2Mx + a 2 . 



(39) 



(40) 



where r m i n (M, a) is the radial tbl value that corresponds to £ = 0. This is chosen as to make $(oo) = ir/2 so that 
both £ and tbl become simultaneously infinite. With this substitution < £ < oo, r nlin < r B L < oo while £(?*bl) is a 
purely monotonic function. 

Performing this coordinate transformation on eqn. (|38[) we obtain 



. , , d 2 R _ , , ,, dR 

A(r BL )-- T - + 2(rBL-M)- 



dr B L 



F(r B h)R = 0, 



with 



^(r BL ) = 



a 2 (oj 2 a 2 tan 2 $ + nMa tan $ + m 2 cos 2 $ - E) 
A(r B L) cos 2 $ 



(41) 



(42) 



Eqn. (|4Tj) has the same 'symbolic' structure as the radial scalar Teukolsky equation in Boyer-Lindquist coordinates. 
The residual difference appears in the form of the F(rBh) function between the two cases. For the Kerr scalar wave 
equation this function is [2| 



[( r BL + a K V ~ a Kfn 

A(r BL ) 



2aYjnuJ — a K ui — E, 



(43) 



9 



where ax is the usual Kerr spin parameter. 

The functions F and Fjz are clearly different; this is not surprising given the fundamentally different physical 
character of the Euler and Kerr fields: (i) Black holes have an event horizon which acts as an one-way membrane for 
impinging waves. No such boundary exists in the Newtonian problem 6 , (ii) The Kerr metric contains rotation- induced 
frame-dragging and this has a distinct impact on prograde and retrograde propagating waves. This effect is encoded 
in the linear m terms in Fk ■ The Euler field, being oblivious to the direction of orbital motion, cannot discriminate 
between prograde and retrograde motion, and that is why F(r) features m? terms only. The Kerr-Euler dissimilarity 
can also be exposed by taking the double limit a — > 0, ok — > at which the two fields collapse to the Schwarzschild 
and Kepler fields respectively, which are fundamentally different. 

All these differences between the Kerr and Euler fields, regardless of their nature, refer to strong-gravity properties. 
Hence, there is still hope that F and Fk (and as a result the wave equations themselves) could approach each other 
asymptotically far from the gravitating source. In fact, this should be required since both fields should approach the 
Keplerian potential at tbl — > oo. Indeed, it is not too difficult to verify that in this limit we have F, Fk — > w 2 r BL . 

The next order asymptotic term of F depends explicitly on the scale parameter k. The choice K = Auj 2 leads to the 
identical behaviour 

F, F K « lo 2 tI L + 2Mtu 2 r BL + 0(r BL ) (44) 

at large tbl- However, this is the furthest one can push the match between F and Fk- Once higher order terms are 
introduced, the two functions begin to deviate. 

To conclude this Section, we return to the two basic Euler wave equations (|37p and (|38p and briefly discuss 
their solutions. As we mentioned earlier, the former equation is solved by the s = spin-weighted spheroidal 
harmonics 2]. The latter equation is solved by the generalized radial spheroidal harmonics of the first and second 
kind EE™*- 1,2 -* (wa, inaM (for a detailed discussion of these functions see [lH). The computation of such special 
functions though, especially in the case of imaginary arguments, is an intricate issue (cf. [HI) and thus it will not be 
attempted here. 

4. CONCLUDING REMARKS 

Analogies between diverse physical systems is a much cherished property of Nature. In this paper we have discussed 
one such analogy between an extremely important, from the astrophysical point of view, relativistic solution of the 
vacuum Einstein equations and a simple Newtonian potential that was first studied more than two centuries ago. 
The Euler field is the unique Newtionian potential that shares the key properties of the Kerr metric, namely, the 
separability of the Hamilton- Jacobi equation (which leads to a Carter-like constant) and the scalar wave equation. 
At the same time the Euler and Kerr fields share the same recurrence relation for their respective mass multipole 
moments. 

As far as the wave equation in the Euler field is concerned we have shown that the equations produced after the 
separation of variables can be cast in a Kerr-like form. While the angular equations are in fact identical the same is 
not true for the radial equations. The two equations coverge to each other only in the weak gravity asymptotic region 
of large radii. Closer to the gravitating source (strong gravity region) the deviation is significant, an expected result 
given the absence of event horizons and frame-dragging in Newtonian gravity. 

Future work on the subject will attempt exploiting the striking Euler-Kerr analogy at a more practical level, for 
instance, by studying the orbital motion of particles in the Euler field, and using this analogy to infer orbital properties 
of realistic cases related to a Kerr black hole as a central object. 
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It can be shown that r m ; n (M, a) is located outside the corresponding Kerr event horizon, r m ; n > M + \/ M 2 — w 
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Appendix A: The solution of V 2 Ki Cp = 

In this appendix we show that the only separable gravitational field in spheroidal (either prolate or oblate) coor- 
dinates is the Euler two-center potential. Considering the potential (|16p , we assume the following expansions for the 
free functions F(£),G(rf): 

+00 +00 
^(O-E^""' G(r,) = J2PkV k (Al) 

k=-l fc=0 

These power-series ensure that the two functions are well-behaved at spatial infinity £ — ¥ 00 and on the equatorial 
plane 7/ = 0. Then 

P 1 + °° 

v «> = + 7« E (»r* + ) ■ (A2) 

^ 1 fc=0 

The first term in this expression solves the Poisson equation identically (the Laplacian operator is given in eqn. (|2ip). 
The remaining terms lead to 



+00 



v 2 y scp = -> E 



k=0 



(k - l) PkV k - 2 { (k - 2)r, 4 + ry 2 [4 — fc =F (2 + fc)£ 2 ] ± ki 2 } 



(k + l)q k C {k+2) { (k + 2)£ 4 T ^ 2 [k + 4+(k- 2)ry 2 ] + krj 2 } 



0. (A3) 



This series vanishes only if q k = for all k > and pk = for all k > except of k — 1. As a result, the full solution 
takes the desired Euler form 

where are constants. 

Appendix B: Multipole moments in spheroidal coordinates 

In order to calculate the specific mass moments of Ve (or V b e) in spheroidal coordinates we first write 

Mf °° 

Ve = U2 2T = y2 a 2nQ2n(0 P 2n(v), 

^ sE = 7e2^ 2\ = y2 a 2niQ2n(iO P 2n(v)- ( B1 ) 

The objective now is to calculate the corresponding coefficients a 2n . To this end we take the product of these 
expressions with P 2 i{rf) and integrate: 

M f +1 P 2l ( V ) 2a 2l 
4 / aj]— - = —— -Qu{Q, (B2) 

M f +1 P 2l ( V ) 2a 2l 

With the help of the integral representation of Qk (Eq. 8.3.3 of [iH ]) 

Qk(z) = \ ( + dt^-, (for Re(z) > -1, lm(z) = 0), (B4) 
Z J _i z — t 
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and the reflection property Qk{z) = (— l) k+1 Qk(— z), the left-hand sides of (|B2[) and (IB3|) can be written as 

2M 



M 
2a 



M 
2ia 



1 C + 7-i f-r? 

and 



drj- 



drj- 



_i ^ + ?7 J_i i£-r? 
Combining these with (IB2[) and (|B3[) leads to 

«2i = -(4/ + l)M/a 
which, remarkably, has the same form for both oblate and prolate potentials 



-0a(0, 



2M 
ia 



(B5) 



(B6) 
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